Ferroelectric Phase Transitions in Films with Depletion Charge 
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We consider ferroelectric phase transitions in both short-circuited and biased ferroelectric- 
semiconductor films with a space (depletion) charge which leads to some unusual behavior. It 
is shown that in the presence of the charge the polarization separates into 'switchable' and 'non- 
switchable' parts. The electric field, appearing due to the space charge, does not wash out the phase 
transition, which remains second order but takes place at somewhat reduced temperature. At the 
same time, it leads to a suppression of the ferroelectricity in a near-electrode layer. This conclusion 
is valid for materials with both second and first order phase transitions in pure bulk samples. In- 
fluence of the depletion charge on thermodynamic coercive field reduces mainly to the lowering of 
the phase transition temperature, and its effect is negligible. The depletion charge can, however, 
facilitate an appearance of the domain structure which would be detrimental for device performance 
(fatigue). We discuss some issues of conceptual character, which are generally known but were 
overlooked in previous works. The present results have general implications for small systems with 
depletion charge. 

77.80.-e, 77.80.Bh, 77.80.Dj, 85.50.-fk 



I. INTRODUCTION 

Ferroelectric perovskite materials used in memory ap- 
plications usually behave as semiconductors It is 
expected, therefore, that at the metal-ferroelectric con- 
tact a depletion zone and the corresponding (depletion) 
charge distribution is formed, similar to situation with 
usual semiconductor-metal contacts. The corresponding 
electronic structure and related effects were extensively 
discussed for layered metal-semiconducting ferroelectric 
systems |^-|^ . The depletion effects are especially impor- 
tant in switching devices . The semiconducting behav- 
ior of some common ferroelcctrics (PZT) was exploited in 
a series of papers [p[-pT| to account for some specific fea- 
tures of switching, fatigue, and rejuvenation phenomena 
in ferroelectric thin films. 

The idea of possible importance of the depletion charge 
for the properties of thin ferroelectric films seems fairly 
natural. However, as will be shown below, ferroelectric 
properties in those systems are very peculiar. The con- 
sistent treatment demonstrates, in particular, that the 
effect of the depletion charge on thermodynamic coer- 
cive field is minute. This is a consequence of long-range 
Coulomb (depolarizing) field which appears simultane- 
ously with the inhomogeneous polarization. The depo- 
larizing field makes the problem of finding the polariza- 
tion field and the point of phase transition non-local, 
and the local, and to a considerable extent, global di- 
electric response rigid. The effects of the depolarizing 



field were apparently neglected by the previous authors, 
Refs. ]9|-pT|]. Therefore, their general relevance is ques- 
tionable because, as is shown below, the depolarizing 
field plays crucial role in determining the dielectric re- 
sponse and the very character of the phase transition. 
The depolarizing field leads to some important phenom- 
ena, like lowering of the critical temperature, Tc, possible 
monodomain-polydomain transition, etc. 

In the present paper we describe the properties of fer- 
roelectric films with the depletion charge in terms of the 
phenomenological Landau-Ginzburg-Devonshire (LGD) 
theory of ferroelectricity (see, e.g. [Q). We restrict our- 
selves to the case of a ferroelectric film with a polar axis 
perpendicular to the film plane. It is assumed that the 
charge distribution remains the same both in the para- 
electric and the ferroelectric phases. 

The key point in our treatment is that the polarization 
is divided into built-in (non-switchable) and ferroelec- 
tric (switchable) parts. Within the paper this division 
is mathematically strict: the built-in polarization, Pbi^), 
corresponds to an extremum of the LGD thermodynamic 
potential. The extremum corresponds to a minimum in 
the paraelectric and a maximum in the ferroelectric phase 
(the latter being similar to a zero polarization solution 
for a pure short-circuited crystal). We shall call the fer- 
roelectric polarization the difference 

Pf=P- Pb- (1) 
The ferroelectric polarization is switchable: there always 
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are two equivalent solutions for -P/(r) in the ferroelec- 
tric phase. In the case of a non-zero bias we include 
the change of the polarization due to the bias into P/, 
whereas the built-in polarization remains unchanged. 

It is convenient to start with a short-circuited film of a 
linear dielectric with a homogeneous density of the space 
charge, p. The electric displacement D in the film is given 
by the electrostatic equation divl? = p, or simply 

dD/dz — p — const, (2) 

where z is the coordinate perpendicular to the film plane. 
Taking into account that the voltage across the short- 
circuited film is zero, one finds: 

D={z-d/2)p. (3) 

From 

D = eoE + P, (4) 

and 

D = eeoE, (5) 
one obtains for the polarization P: 

P=^—^D. (6) 

£ 

Now we shall consider a ferroelectric in a paraelectric 
state. We shall neglect, for a moment, the nonlineari- 
ties and make some qualitative remarks. It follows from 
Eq. @ that P k, D when £ > 1. Note that the electric 
displacement is always determined by the charge den- 
sity, it does not depend on the dielectric constant [cf. 

]. When the phase transition is approached, the di- 
electric constant increases, but the polarization remains 
practically unchanged, unlike the electric field which ex- 
hibits a strong temperature dependence. At the Curie 
temperature (e = oo) the field is zero in the linear ap- 
proximation, while the polarization is almost the same as 
it was far from the transition. In other words it is more 
appropriate to talk about a built-in polarization rather 
than a built-in electric field while considering the effects 
of space charge in ferroelectrics. This conclusion is not 
limited, of course, to the case of homogeneous distribu- 
tion of the space charge. For a medium with a large 
dielectric constant the value of the polarization is close 
to the value of the displacement field (^). Hence, near 
complete compensation of the space charge (produced by 
ionized donors, etc.) by the bound charges, correspond- 
ing to inhomogeneity of the polarization, occurs in a 'soft' 
dielectric medium. 

The difference between the ferroelectric and non- 
ferroelectric (built-in) polarization can be clarified by 
considering effect of the bias voltage. One finds for a 
linear dielectric 



D={z- d/2)p + eeoU/d, (7) 

where U is the bias voltage and the polarization P natu- 
rally splits into two parts, the switchable ('ferroelectric') 
polarization Pf and non-switchable ('built-in') polariza- 
tion Pf,: 

P = Pf + Pt; (8) 
Pf = {e- l)eoU/d, (9) 

P,{z) = {z-d/2)p'-^. (10) 

£ 

Thus, in the presence of both the space charge and the 
bias voltage, there are two contributions to the polar- 
ization, but only one of them 'feels' the phase transi- 
tion. The built-in polarization does not change with the 
bias, it remains inhomogeneous and almost insensitive 
to the phase transition [since (e — l)/e « 1]. The in- 
duced (ferroelectric) polarization Pf — {e ~ 1)£qU /d = 
(£ — l)£o-Ecxt is homogeneous. The corresponding sus- 
ceptibility X/ = dPf /dEcxt o: e — 1 diverges at the tran- 
sition temperature in the linear approximation. We pro- 
vide more details and discussion about the division of the 
polarization into ferroelectric and non-ferroelectric parts 
below. 

The paper is organized as follows. In Sec. II we analyze 
the problem of phase transition in a short-circuited film 
with homogeneous space charge. We assume that ferro- 
electric material is having a second-order phase transi- 
tion in the bulk. We shall show that the phase transi- 
tion in the film with depletion charge is inhomogeneous, 
with both 'global' and 'local' phase transition tempera- 
tures being lower in the film than in the bulk. Thus, the 
presence of the space charge will be shown to somewhat 
suppress the ferroelectricity. 

To treat the ferroelectric phase, in Sec. Ill we consider 
a special exactly solvable case of the space charge distri- 
bution, when all the charges are located at the central 
plane of the film. We find that the presence of the deple- 
tion charge leads to certain lowering of the critical tem- 
perature. The corresponding lowering of the thermody- 
namic coercive field is shown to be minute. Therefore, the 
possibility of a large effect of the built-in field on switch- 
ing in monodomain samples, asserted in Refs. |p|-|rH, is 
rather questionable, especially since they have not con- 
sidered the actual nucleation processes. Other types of 
the charge distribution are discussed qualitatively. In 
Sec. IV we show that in materials with first order phase 
transitions the depiction charge also effectively decreases 
the dielectric constant of paraelectric phase. We summa- 
rize the results in Sec. V and add there some comments 
of conceptual character. 

II. PHASE TRANSITION IN THE PRESENCE OF 
HOMOGENEOUS SPACE CHARGE 
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A. Main equations 

The nonlinear electric equation of state has a usual 
Landau form: 



AP + BP^ = E, 



(11) 



where P{E) are the z-components of the polarization 
(electric field). We assume, as usual: 



a{T -Tc), B = const, 



(12) 



where Tc is the Curie temperature and a =:const. To- 
gether with Eqs.(0),(0) and the boundary condition 



Edz = [/, 



(13) 



we now have the complete set of equations to treat the 
ferroelectric phase transition in an electroded film. Note 
that we cannot use Eq. (|^) now because it implies a linear 
relation between D and E, whereas we are interested in 
a proper account for the non-linear effects. 

Integrating Eq. (0), we find with the use of M): 



£qE{z) + P{z) = pz + K, 



(14) 



where K is & constant which should be found from 
Eq. (|l|). We obtain: 

eQE{z) = -P{z) -f (z - d/2)p + i y P{z)dz + eo-Bext, 

(15) 

where the external field is produced by external bias volt- 
age, 

E,y,t = U/d. 

Combining this with Eq. ( |Tl| ) we get: 

(A + £oi)P + Bp3 _ 1 j p(^)rf^ (16) 
d£Q Jq 

- (z-d/2)p/eo + ^cxt. (17) 
From Eqs. (|l^),(|l|) we obtain, up to terms oc P^ : 



{A + 3BP^ + eo^)Pf + 3BPbPf + BPf 



= E,, 



(18) 



B. Paraelectric phase 

Let us first calculate P/ in the paraelectric phase to 
the linear approximation, i.e. for small U . We have to 
solve the simple integral equation 



{A + 3PP2 + £^^)Pf --^ [ Pf{z)dz = Sext, (19) 
MEO Jo 



assuming that A is positive. Introducing a notation 



we find from (|l9|): 



P, 



EcKt + I 



f 



A + 3BP^+eo^' 



(20) 



(21) 



To find / we substitute this expression into Eq. ( P0[ ) and, 
as a result, find P/ in an explicit form: 
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Pf{z) = E,,,^—^ 



1 



A + 3BP^{z) + So 



where 



Thus: 



G = 



1 

dso 



dz 

A + WP^ 



(22) 



(23) 



1-G = 



deoJ^ A + 3BP^+eo 



2 ^ ^^jdz^ A + 3B{Pk',, 



(24) 



where {...)= d ^ dz . . . , and we have taken into ac- 
count that typically s^^ > {A,3BP^). 

We see, therefore, that the susceptibility Xf = 
dPf / dEcxt oc (1 — G) ^ diverges at the temperature Tc, 
which is lower that the Curie- Weiss temperature Tc in 
pure crystal (Fig. |l|), and the divergence of the dielectric 
response would take place when A = —3B{P^), i.e. at 
T = fc: 



Tc-Tc 



3B{P^) 



(25) 



Let us estimate this difference. For a displacive transi- 
tion B and a have normal 'atomic' values ||l2|. For P^ w 
Pat, where Pat is the atomic polarization, Tc—Tc would be 
about Pat ^ 10'*— lO^K. For the charged impurity concen- 
tration Nd one estimates Ph ^ {Nd/Nat){d/ds,t)Pat where 
iVat is the 'atomic' concentration (10^^ — 10^^)cm~'^, dat 
is the 'atomic' distance (unit cell length, 3-5A) With the 
donor concentration Nd — lO^^cm"'^ the thickness 
of the totally depleted film (equal the width of the de- 
pletion layer) would be about 2-10'^A, and we estimate 
Pb - (10-2 - 10-3)p^t and Tc-fc {I - 100)K. We 
see that the phase transition temperature can be reduced 
substantially, i.e. the paraelectric phase can be consider- 
ably 'harder' than the same phase in the pure crystal. 

We see already that there is a radical difference be- 
tween the effects of external field (the bias voltage) and 
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perature from Tc in pure sample down to Tc by (1 — 100)K 
(see text). 



the built-in field (polarization) on the ferroelectric tran- 
sition. The former would smear out the phase transi- 
tion and shift the maximum of the dielectric response to 
higher temperatures, while the built-in field does not lead 
to the maximum in the dielectric response at A > 0. 



FIG. 2 Bratkovsky and Levanyuk 

FIG. 2. Schematic of the dielectric response 

Xj^ = dEcxt/dP which becomes inhomogeneous in the pres- 
ence of the depletion charge (a), and the built-in polarization 
Pt (b). 



C. Ferroelectric phase 

But what happens at A < 0? This is not so easy to 
determine. We shall try to answer this question, at least 
semi-quantitatively, but first we have to explain what 
the difficulty is. It is seen from Eq. ( p2|) that Pf is inho- 
mogeneous over the sample, being smaller near the elec- 
trodes. It is worth mentioning that inhomogeneity of 
Pf under the bias voltage is different from the inhomo- 
geneity of the dielectric susceptibility of the paraelectric 
phase. Indeed, the 'local' dielectric susceptibility is what 
one would find by cutting out a small piece of the fer- 
roelectric at a given position in the sample, making a 
capacitor out of this piece, and measuring its dielectric 
response. In our case we have to 'fill up' the capacitor 
with a medium with a given local values of Pb'. i.e. it is 
the built-in polarization that makes the ferroelectric re- 
sponse of the film effectively inhomogeneous. If P^ were 
constant, the Eq. (|l^) would have had a homogeneous 
solution: 



A 



>> A in the latter case. The smaller 



{A + iBP^)Pf^E, 



(26) 



One can see that the local susceptibility xi 



iBPi) 



dPf/dE,,t = {A 
than the polarization Eq. ( 
3BP^{z) in the case of xio 



is more inhomogeneous 
): one has to compare 
where A is small, with 



inhomogeneity of the polarization Pf in relatively more 
inhomogeneous dielectric medium (xioc) is due to the mu- 
tual Coulomb interaction of the bound charges ('longitu- 
dinal' inhomogeneity of the polarization). Consequently, 
the system tries to lower this energy by reducing the in- 
homogeneity.Therefore, the phase transition takes place 
in a fairly strongly inhomogeneous medium. There is 
also another inhomogeneity, which is described by the 
second term on the left hand side of the Eq. (18). It is 
expected to be important in the ferroelectric phase. 

In the center of the film, where Pb = 0, the local sus- 
ceptibility is the same as in pure crystal but at other 
points of the crystal it is 'harder' (meaning that xioc is 
lower), the hardest region being close to the electrodes 
(Fig. H). If not for long-range Coulomb interactions, this 
would means that the ferroelectric phase transition would 
occur at the critical temperature for pure sample, T^. 
However, the Coulomb interaction suppresses the phase 
transition across the system, ans Tc reduces down to Tc, 
Fig. 1. 

Since the problem of the film with depletion charge 
proves to be fairly complicated, it is instructive to con- 
sider first some special cases, where the treatment is eas- 
ier, and get a feel of the relevant effects. The most natural 
way to simplify the problem is to replace the continuous 
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III. PHASE TRANSITION IN THE FILM WITH 
SPECIAL FORMS OF THE SPACE CHARGE 
DISTRIBUTION 

A. Charge located at the central plane and its effect 
on the coercive field 

We denote by Pi, Ei and P2, E2 the polarization and 
the field in the two halves of the short-circuited film with 
the external bias voltage U (Fig. ||). We then have: 



2U 



2E, 



(27) 
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eoE2 + P2- eoEi - = a, 



^2f 



P2, 



FIG. 3. Bratkovsky and Levanyuk 



(b) 



FIG. 3. (a) Schematic of the exactly solvable case when 
the charge (with surface density a) is located on a plane 
in the middle of a ferroelectric film under bias voltage U. 
The built-in polarization Pi, takes the values —{+)Po, where 
Po = o"/2, in the left (right) halves of the film (regions 1 and 
2, respectively), (b) Possible domain structure with ferroelec- 
tric (switchable) polarization which has different magnitude 
in the left and right halves of the film. 



distribution of the built-in polarization by a stepwise dis- 
tribution. The simplest problem of that type would be a 
charge located at a plane in the middle of the film with a 
constant charge density a, Fig. ||. In this case the built- 
in polarization has opposite signs but the same absolute 
value in the two halves of the film, so that the film is 
homogeneously 'hard'. It is not completely homogeneous 
though, there is the inhomogeneity due to the second 
term on the left hand side of Eq. (^8|) . But now we can 
study the effects of two inhomogeneities separately, be- 
ginning with the case of a 'harmonic' homogeneity and 
a weaker 'anharmonic' inhomogeneity. We would expect 
that in this case the phase transition occurs into a mon- 
odomain state and show below that this is a justifiable 
assumption. 



(28) 



APi(2)+SPi%) =i?i(2)- (29) 

Finding Ei, E2 from the first two equations, one can 
present Eq. (p9|) in the form: 



1 



,3 ^2 ^ 



2eo / 2£o 2eo 



A 



1 

2^ 



P2 



BP^ - — = Eo 
' 2en 



2^' 



By summing up these two equations we get: 
(Pi + P2) [A + B{P^ - P1P2 + P|)] = 2E,; 



(31) 



(32) 



Obviously, for the case of unbiased sample (Pcxt = U = 
0) there is always a solution Pi = — P2 = — Pq (Fig. ||) 
where Pq is given by 



-en')Po 



BP^ = a/2eo. 



(33) 



Since Pq is much smaller than the saturation polariza- 
tion Ps , and accounting for the fact that in ferroelectrics 
^1 <C (it reads as \A\ <C 47r in CGS units) the solu- 
tion is 



(34) 



Indeed, the estimate for a typical PZT film (see below) 
gives Pq — 2.5/iCcm^^, whereas the saturation polariza- 
tion Pg = 40/iCcm^^. We shall consider this solution as 
the built-in polarization and present the other solutions 
in the form: 



Pi{2) — -{+)Pa + Pi{2)f, 



(35) 



where Pi{2)f is the 'ferroelectric' polarization. From 
Eqs. (10), (ml) we have: 
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A + 3BP^ 



2er 



1/ 



2eo 



A + 



1 



2/ 



2eo 



2£o 

12 I D d3 



3BPoP^f + BP^f ^ E,^f 



(36) 



(37) 



These equations can be considered as obtained by mini- 
mization of the thermodynamic potential: 



1 

2F0' 



PfP 



B 



-BPo (Pf^ - P|^) + -(P4^ + P4^) 



-i^cxt (P1/ + P2/), 
where we have introduced 

A = A + WPo. 



(38) 



(39) 



It is instructive to go over to new variables, the average 
switchable polarization and the difference between the 
switchable polarizations in both halves of the film: 



P 



_ P1/+P2/ 



AP = Pi/-P2/. 



(40) 



Now we are in a position to find the ferroelectric polar- 
ization in the ferroelectric phase. The equilibrium values 
P, AP at A < are: 



p2 ^ p2 ^ 

AP = 6soBPaP^ 



-6eoAPo. 



(42) 



These formulas are approximate: we have neglected an 
inessential renornialization of the coefficient B due to 
the coupling given by the fifth term on the right hand 
side in Eq. (|l|), and we have taken into account that 
I A \« Eq^. Then the solution for spontaneous ferro- 
electric polarization is: 



Pi 



1/ 



P 



2/ 




±P, 



(43) 



3eo\A\Po = ±P, - -eo\A\a, (44) 



where -l-(-) sign corresponds to right- (left-) ward directed 
polarization in the domains, Fig. 3(b). We see that in the 
half of the film where the ferroelectric polarization has 
the same direction as the built-in polarization (right half 
in Fig. H), its value, P2/, is smaller than in another half. 



In these variables the thermodynamic potential takes the 
form: 



F = AP^ 



B 



2 

A + Eo' 



-P" - 2ScxtP 



Ap2 - SBPnP^AP 



■^BPqAP'^ + — AP* + -BP^AP^. 



(41) 



The stability of the paraelectric phase (Pi/ = P2/ =0) 
will be lost a,t A = A + 3BPq = 0, therefore, the param- 
eter P is 'critical' at the phase transition while AP is 
'non-critical', since the coefficient before AP^ docs not 
go to zero. Now it is evident that the phase transition 
we discuss is second order, since there is no cubic term 
in P. The external bias voltage couples only to P, so 
that the phase transition for the order parameter P will 
be smeared out by external field, whereas it has no effect 
on the built-in polarization Pq. This highlights the qual- 
itatively different response of the system to external and 
built-in field due to space charge. 

We see from the third and fourth terms in F that 
the equilibrium value of AP is proportional to P^ and, 
therefore, the last three terms in (41) would implicitly 
contain terms higher order in P, like P^. This would be 
inconsistent with the form of the Landau functional that 
we used initially. Therefore, the last three terms in ( pl| ) 
should be omitted. 



B. The effect of the depletion charge on a coercive 
field 

Now we are able to estimate the effect of the deple- 
tion charge on a coercive field, and it happens to be very 
small. From the free energy (^ij) we see that the depen- 
dence of the average switchable polarization P on i?oxt 
has a form of a hysteresis loop 



AP + BP^ 



(45) 



Here we accounted for the fact that the fifth term in the 
free energy (^ is actually oc P'^ because of (^, and 
this is equivalent to the renormalization of the coefficient 
B ^ B = B{1 - ISeoPPo). 

The 'thermodynamic' coercive field for the hysteresis 

loop P = P(Sext) is 



Er = 



-A 



3/2 



33/251/2 



(46) 



One can easily find the hysteresis loops for Pi / and P2 / 
from Eqs. (0),(|2|): 



Pi/ = He. 

P2/ = P{E. 



cxt J 
cxt ) 



3eoPPoP'(Scxt) 
3eoPPoP'(Scxt) 



(47) 



The hysteresis loops are interesting since the thermody- 
namic coercive field is exactly the same in both halves 
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donor concentration Nd = 10"'^*cm 




FIG. 4. Schematic of the hysteresis loop P = P(-Ecxt) of 
the film with the charge at the central plane, Fig. 3. Note 
that the thermodynamic coercive field Ec is exactly the same 
in both halves of the film (as seen from the hysteresis loops for 
Pif and P2/). The polarization is normalized to Pc = P{Ec). 



ously i. 
Po = <j/2 



, discussed previ- 
With these parameters a = NdL = 5^C/cni^, 
= 2.5/LtC/cm2. We then find from (||) 



= 0.02 



Obviously, the effect quickly gets even smaller for thinner 
films of practical interest, lS.Ec cr^ (x . 

Thus, we see that the effect of the depletion charge 
on the coercive field exists, but it is negligible. It has 
been speculated that the 'experimental' coercive field for 
nucleation Ecn should be reduced by exactly the value of 
the built-in field E},, Ecn = £^oxt + E^ Since these 

authors have estimated E), « lOOkV/cm, whereas the 
observed coercive field was Ec ~ 80kV/cm |], one would 
conclude that the coercive field has been suppressed by 
more than a half by the built-in field. However, this result 
and the whole previous analysis are rather questionable, 
since no actual nucleation processes have been considered 
in Refs. H-jlll. 



C. On the possible domain structure 



of the film. Fig. in spite of that difference in absolute 
values of the polarizations. 

Since in the ferroelectric film with depletion charge 
\A\ < \A\, the coercive field there is smaller than the one, 
Eco, in a sample without the space charge. The physical 
reason for this reduction is the lowering of the transition 
temperature at the phase transition point, where the co- 
ercive field is zero. Thus, one can expect that the switch- 
ing would be somewhat easier in the material with the de- 
pletion charge in comparison with a pure material. This 
observation is only suggestive since the actual switching 
occurs at much lower fields and proceeds by nucleation 
and growth processes, whereas the 'thermodynamic' co- 
ercive field Ec refers to thermodynamic limits of stability 
of the homogeneous polarization in the external field. 

Let us estimate the change in the coercive field for 
typical thin PZT film not very close to the transition, 
|A| > 3BP^, from (H) 



where 



E. 



cO 



AEc 



2\A 



Ec 



3/2 



Eco — AEc 



(48) 



33/251/2^ 



-^Ec,^{l-2eo\A\) 



q p2 

-^^co^, (49) 



where Eco is the coercive field in the pure material, 
and we again used the fact that |^|eo <C 1. We now 
consider the PZT film of the thickness L — 300 nm, 
the saturation polarization Pg = 40 fiC/cm^, and the 



One sees from Eq. (^) that there is a discontinuity 
of the ferroelectric polarization at the central plane. It 
creates depolarizing field, which can, in principle, be re- 
duced by formation of a domain structure [Fig. §(b)]. To 
see if it really takes place we have to use the LCD ther- 
modynamic potential with the energy of the electric field 
included explicitly: 



LGD 



= / dV 



A 



B 



-P" 



K (dP\ EqE^ 
2 \'dx ^ 



(50) 

where the gradient term is added because we would like 
to consider a possibility of domain structure formation 
with the domain walls perpendicular to x direction, and 
this term gives rise to a surface energy of the walls. 
Once more we divide the polarization into two parts as 
in Eq. (m) and write 



E = Eb + Ef, 
with obvious notations. Then 



LGD 



Fi + Fy+ F-i 



(51) 



(52) 



where 



Fi 



dV 



A + 3BP^ , o B , 

2 



^2\dx J ^ 2 



(53) 



7 



= J dViEbPf+eoEbEf)^ J dVEbDf, (54) 



F3^F3{Pb,Eb 



(55) 



The energy F2 = for the problem of calculating the 
energy of the domain structure periodic along the x axis, 
and we can consider Fi only. The domain structure is 
shown in Fig. ^. We suppose that the domain width is 
a < d, otherwise the domain structure would not form. 
The electric field is created by the bound charges at the 
central plane of the film. This field is approximately 



1 - - '^A 



(56) 



and it is concentrated in a layer of width « a. The energy 
of this field per unit area and unit length along x direction 
is 



eoA^Pga, 



(57) 



where we have omitted a numerical factor. The surface 
energy of the domain walls for the same region is 



jd 



(58) 



where 7 is the surface energy of the domain wall. It is 
known to be 



7 



B ^' 



(59) 



where ^ is the domain wall width. Minimizing the to- 
tal energy, one finds the period of the suspected domain 
structure: 



d^ 



BPS 



1/2 



The condition a < d is fulfilled if 



d> 



P. 



BPg 



Pr 



i2 ^ 



(10^ - io^)e 



(60) 



(61) 



We have already mentioned that Pat/Po ^ 10^ — 10'^, so 
that the prefactor of is about 10^ — 10^. Taking into ac- 
count that f is no less than the unit cell length, and that 
it can easily be an order of magnitude larger than that, it 
seems doubtful that the formation of the domain struc- 
ture due to discontinuity of the polarization in the center 
of the film would be possible in films with d ^ lO^dat. 
In any case, if it forms, then at fairly low temperatures, 
so that our initial assumption that the phase transition 
takes place into a monodomain state is justified. We see 
that the 'anharmonic' inhomogeneity is not the most im- 
portant reason for the domain structure formation, and 
we shall neglect it in the future. 
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FIG. 5. Bratkovsky and Levanyuk 



FIG. 5. Schematic of a ferroelectric film with two charged 
planes inside. The charged planes lead to piece- wise built-in 
polarization Pi,. 



D. Charge located at two symmetric planes 

A stepwise distribution of the built-in polarization, 
which is more like that in Fig. ^, can be obtained with two 
similar symmetrically positioned charged planes (Fig. ^). 
We have now two 'hard' near electrode regions and the 
central 'soft' one. For A + '^BP^ > the system re- 
sembles a plate of a pure ferroelectric with two 'passive' 
layers close to the electrodes. Phase transitions in a sim- 
ilar system (ferroelectric plate with vacuum gaps near 
the electrodes) were studied by Chensky and Tarasenko 
[ p^ . They found that there are transitions to both mon- 
odomain and polydomain states depending on geometri- 
cal parameters. Generalizing their consideration to the 
case of two dielectric layers we found that the transition 
to the monodomain state takes place if 



di < die = £1 



12£oK 



1/2 



(62) 



where di, ei are the width and the dielectric constant 
of the 'passive' layer, e± is the dielectric constant in the 
plane of the film (perpendicular to easy axis) . In our case 
the value of die d can be estimated as: 



di 



p2 
^0 



(10^ - 10^)1. 



(63) 



This estimate is certainly crude but it seems likely that 
the phase transition in completely depleted films pro- 
ceeds into a monodomain state. We have recently dis- 
cussed a domain structure in a ferroelectric plate with 
'passive' layers |jl^. It was found that far from the phase 
transition a domain structure exists for however thin pas- 
sive layer. Therefore, one can expect that a domain struc- 
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ture forms at some temperature below the phase transi- 
tion in a short-circuited completely depleted ferroelectric 
film, but this question should to be studied in greater de- 
tail. 



IV. FIRST ORDER PHASE TRANSITION 

It might be tempting to conclude, based on Eq. (p5|), 
that while the built-in polarization 'hardens' ferro- 
electrics with second order transition it 'softens' those 
with first order transitions, since the coefficient at in 
the LGD thermodynamic potential is negative in the lat- 
ter case {B < 0). However, such a hastily conclusion is 
almost certainly wrong. 

The fact of the matter is that, while considering rel- 
atively 'weak' first order phase transition (otherwise the 
LGD theory is not applicable), one has to take into ac- 
count explicitly the interaction of the order parameter 
with elastic deformations. Specific features of the elas- 
ticity in solids make the above statement about the sign 
of the coefficient B before P'^ fairly vague. There are 
several different coefficients and some of them remain 
positive for weakly first order transitions jl^]. To treat 
the problem, one has to start with the LGD potential in 
the form: 



A 



B 



C 



F = -P^ + —P" + —P" 



6 



K 



+rP uii + —Uii + Uik ~ -^^ikUii 



1 



(64) 



where Uik is the strain tensor, r is the electrostriction 
constant, K , jjL are the bulk and the shear moduli. For 
a bulk sample the phase transition takes place usually 
in a free crystal and after minimization over the strain 
components one obtains: 



F=ip^ 



B 



.p4 



C 



(65) 



tensor is u^z- Putting all other strain components in 



Eq. 



where 



to zero and minimizing over Uzz we obtain: 



F^^P'^ + ^P*. 
2 4 



C 



B = B 



2r2 



K 



1^ 



^ B 



2r' 



K K 



(67) 



(68) 



It is worth mentioning that K + \s the modulus with 
respect to uniaxial extension without changing the lateral 
size 10 . The last factor in Eq. (|6^) is usually of the order 
of unity and if | | is not very large then _B > and the 
results of Sec. II are applicable for the case considered 
there as well. 

One might conclude that, because i? > 0, the first or- 
der phase transition in a bulk sample becomes second or- 
der in the film. Such a possibility is not excluded but sev- 
eral reservations are in order. Firstly, we have assumed 
that only one component of the strain tensor operates. 
It is a natural assumption for an infinite film but the real 
films are not infinite. Secondly, even for an infinite film 
there are various possibilities for the phase transition: it 
could, for instance, transform into a mixture of homoge- 
neous and inhomogeneous (two phase) state. Which of 
these possibilities is realized in practice depends on the 
conditions that we have not specified here, e.g. on the 
elastic properties of the substrate. The nature of first or- 
der phase transitions in films is worth studying together 
with taking into account the actual experimental condi- 
tions. 

In conclusion of this Section we would simply note that 
in materials with first order ferroelectric transition the 
paraelectric phase in thin depleted films becomes more 
'ferroelectrically rigid', its dielectric constant decreases, 
in analogy with ferroelectrics with second order phase 
transition. 



where 



B = B- 



1<~ 



(66) 



The first order of the phase transition means B < 0, but 
the coefficient B itself may be positive and it is indeed 
positive and large for not too 'strong' first order struc- 
tural transitions, because 2r'^ / K is of 'atomic' value in 
this case, and B is rather small. 

But for a phase transition in an inhomogeneous sys- 
tem the above procedure does not work. One cannot 
minimize over the strain components because they are 
not independent there |p^ . For our case, one can assume 
that all the displacements are perpendicular to the film 
plane so that the only relevant component of the strain 



V. COMMENTS 

We have found that the polarization in the present 
problem naturally separates into 'switchablc' and 'non- 
switchable' parts in the presence of the built-in (space) 
charge. Obviously, only the switchablc polarization may 
be used as an order parameter in the problem. In our 
example of a system with all the charge placed at the 
central plane the order parameter was P — Po(r) averaged 
over the sample, where Po(?0 is the non-switchable part 
of the polarization. The case of a homogeneous charge 
in the ferroelectric proved to be too complicated and its 
consideration was omitted. What is the problem with 
selecting the order parameter in the present situation? 
We would like to clarify the question in this section. 
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One difficulty with selecting the order parameter in 
systems with depletion charge is that these systems are 
inhomogeneous with respect to the phase transition, even 
when the depletion charge is homogeneous. The order 
parameter for a phase transition in the inhomogeneous 
medium is the amplitude of a function that has a space 
distribution of the quantity that plays a role of the or- 
der parameter in homogeneous medium. The function 
gives the form of the 'most rapidly growing fluctuation' 
or of the symmetry perturbation responsible for the loss 
of stability of a symmetric phase. 

Another difficulty stems from the polarization being 
an electric quantity, intimately related to a distribution 
of the electric charge in a crystal. In defining the order 
parameter one has to recall that in theory the notions 
of the 'order parameter', 'phase transition', and so forth, 
implicitly refer to an infinite medium. It is ambiguous, or 
even meaningless, however, to consider the polarization 
of an infinite medium (see e.g. Ashcroft and Mermin's 
book, Ch. 27). This is because the electric field in a sam- 
ple is determined by the boundary conditions (shape of 
the sample, presence or absence of electrodes, bias volt- 
age, etc.). One can avoid specifying the boundary condi- 
tions and still remain within rigorous definitions by con- 
sidering the plane waves Pk of the order parameter in an 
infinite medium and then taking a long wave length limit 
k ^ 0. Those waves which stiffness constant goes to zero 
at the phase transition (for almost infinite wave lengths) 
belong to the order parameter, otherwise they do not. In 
most cases this procedure is not needed, since there is no 
singularity (discontinuity) at k = 0, so that the values of 
all quantities at k = and at k ^ are the same. This 
is not valid in the case of polarization. Only transversal 
polarization waves can loose stiffness at a critical point 
(i.e. the transversal mode softens to zero at Tc), whereas 
the longitudinal waves, unlike the transversal ones, cre- 
ate a macroscopic electric field and, therefore, remain 
hard at the phase transition. One can identify the or- 
der parameter of a (proper) ferroelectric transition with 
the transversal part of the polarization, but this is not 
convenient in practice. One usually considers finite size 
samples with well defined boundary conditions, and it 
might be quite difficult (if not impossible) to discrimi- 
nate between transversal and longitudinal polarizations. 
The best we can advise at this moment is to remember 
that there may be a part of polarization that does not 
correspond to the ferroelectric order parameter. Direct 
solution of the 'electric' equation of state ( pi] ) together 
with Maxwell equations (|^) with account to boundary 
conditions will do the job. The simultaneous solution of 
the equation of state and electrostatic equations makes 
it clear that, for instance, the effects of the external field 
and the space charge are qualitatively different. 

Usually the complication about selecting the order pa- 
rameter does not arise, but in the present case it is conve- 
nient to discriminate between the ferroelectric and non- 



ferroelectric part of polarization: what we call the 'built- 
in polarization' is clearly 'non- ferroelectric', it is longitu- 
dinal and practically does not change at the ferroelectric 
phase transition. What we call the 'ferroelectric' polar- 
ization is strictly 'transversal' (in the sense that it creates 
no field or, rather, this field is screened by the charges 
at the electrodes) only when this polarization is homoge- 
neous. It is not exactly our case, so that the conceptual 
basis for such a division of the polarization may easily be 
questioned. Therefore, we would refer to Eqs. (|lj),(|l^) 
and our present comments. This division is very helpful 
for the present problem, with simplifications and com- 
plications resulting from the interaction between the two 
parts of the polarization, since for ferroelectrics the 'elec- 
tric' equation of state is nonlinear. The effects of this 
interaction were focus of the present paper. 

VI. CONCLUSIONS 

We have demonstrated a peculiar nature of the phase 
transition in ferroelectric films with depletion (space) 
charge, and specific properties of the ferroelectric phase. 
The charge leads to appearance of the built-in (frozen) 
polarization which is not sensitive to the ferroelectric 
phase transition, and it suppresses the ferroelectricity in 
the near-electrode regions. This is true of both second- 
and first-order phase transitions. The lowering of the 
critical temperature leads to very small reduction in the 
thermodynamic coercive field. The crux of the mater 
is that the system is strongly affected by long-range 
Coulomb field, accompanying the inhomogeneous polar- 
ization, which makes the local and, to some extent, the 
global dielectric response rigid (meaning the reduced di- 
electric susceptibility). The main effect is that it sup- 
presses the critical temperature across the whole sam- 
ple. The effects of the depolarizing field in this prob- 
lem have been apparently neglected by previous authors 
[PI-pl, who also speculated about built-in field assisted 
switching in ferroelectric films. Since no actual nucle- 
ation processes have been considered in these works, we 
doubt those speculations have any justification. 

The unusual feature of the present situation is that 
the value of the ferroelectric (switchable) part of the po- 
larization is smaller when it is parallel to the built-in 
field, and larger when it is anti-parallel. This may facili- 
tate a splitting of the film into domains at low tempera- 
tures, thus playing a detrimental role in the device per- 
formance (fatigue). This illustrates an important point 
that the switching behavior is defined by both the bulk 
properties of the sample (electronic, ferroelectric, and mi- 
crostructural) and boundary conditions, which all should 
be taken into account. The observed, i.e. hysteretic, be- 
havior of the device may change drastically depending on 
variations of any of these conditions and the ferroelectric 
material itself (like the growth of the 'passive' dielectric 
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layer at the electrode- ferroelectric interface). We expect 
that the behavior discussed in the present paper is generic 
for systems of small size where screening of the impurity 
charge becomes ineffective and depletion effect becomes 
pronounced. 
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